Time scale for the onset of Fickian diffusion in supercooled liquids 
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We propose a quantitative measure of a time scale on which Fickian diffusion sets in for supercooled 
liquids and use Brownian Dynamics computer simulations to determine the temperature dependence 
of this onset time in a Lennard- Jones binary mixture. The time for the onset of Fickian diffusion 
ranges between 6.5 and 31 times the a relaxation time (the a relaxation time is the characteristic 
relaxation time of the incoherent intermediate scattering function). The onset time increases faster 
with decreasing temperature than the a relaxation time. Mean squared displacement at the onset 
time increases with decreasing temperature. 

PACS numbers: 64.70.Pf, 61.20.Lc, 61.43.Fs 



Understanding the origin of the extreme slowing down 
of liquids' dynamics upon approaching the glass transi- 
tion and the nature of the transition itself has been of 
great interest for several decades. A lot of recent activity 
has been stimulated by the recognition that close to the 
transition the liquids' dynamics become not only very 
sluggish but also increasingly heterogeneous 1]. While 
the presence of dynamic heterogeneities is commonly ac- 
cepted, the details of their spatial and temporal struc- 
ture have been only partially established. In particular, 
the question of the lifetime of dynamic heterogeneities 
is quite controversial: whereas two separate experiments 
|5|, Ig found that near the glass transition the lifetime 
is significantly longer than the a relaxation time, other 
experimental studies 0) S found the lifetime to be 
comparable to the a relaxation time. In principle, the 
controversy can be resolved by postulating that the tem- 
perature dependence of the lifetime is stronger than that 
of the a relaxation time 0. However, the physical inter- 
pretation of the new time scale remains unclear. On the 
computational side, there have been a few attempts to 
estimate the lifetime of dynamic heterogeneities. Most 
of them 10, 0, 0, 01 found the lifetime to be com- 
parable to the a relaxation time. To the best of our 
knowledge, the only exception is a recent preprint in 
which it is shown that, in a kinetically constrained spin 
model resembling a fragile glass former, the lifetime is a 
few times longer than the a relaxation time. More im- 
portantly, Ref. 14] found that the lifetime increases with 
decreasing temperature somewhat faster than the a re- 
laxation time. It should be noted that one earlier study 
(lif also found that at a low temperature the lifetime of 
dynamic heterogeneities is a few times longer than the 
a relaxation time. However, a careful study of the tem- 
perature dependence of these two times have not been 
performed. Thus, the question of the existence of a time 
scale longer and increasing faster than the a relaxation 
time remains unresolved. 

The goal of our study is to investigate the temperature 
dependence of a different characteristic time that is re- 
lated to the lifetime of dynamic heterogeneities: the time 



tion of the logarithm of single-particle displacement, 
logio(fr), during time t, P{\og^^{8r)-t) [3, [ll, I£, iS] • 



for the onset of Fickian diffusion. 

In order to define the onset time we use, as an in- 
dicator of Fickian diffusion, the probability distribu- 

le displace 

This distribution is defined in such a way that the in- 
tegral J^^ P{x;t)dx is the fraction of particles whose 
value of \ogiQ{5r) is between xq and Xi. The proba- 
bility distribution P{logiQ{Sr);t) can be obtained from 
the self part of the van Hove correlation function '2di, 
P(logio((5r);t) = ln{lO) Air Sr^GsiSr,t). The probability 
distribution P{logiQ{5r);t) is a convenient indicator of 
Fickian diffusion because if particles move via Fickian 
diffusion then the self part of the van Hove function 
is Gaussian and the shape of the probability distribu- 
tion P{\ogiQ{dr)]t) is independent of time. In particu- 
lar, the height of the peak of this distribution is equal 
to ln(10)^54/7r e^^^^ « 2.13 and deviations from this 
value indicate non-Fickian particle motion. We define 
the time for the onset of Fickian diffusion, Tp, as the 
time at which the peak of P{\ogiQ{Sr);t) is equal to 90% 
of its value for a Gaussian distribution of displacements, 
P{\ogiQ{Srmax)',Tp) « 1.92. Wc will discuss the thresh- 
old value of 90% together with a different indicator of 
Fickian diffusion at the end of this Letter. 

It should be noted that a deviation of the probabil- 
ity distribution P{\ogiQ{Sr);t) from its universal shape 
expected for Fickian diffusion indicates dynamic hetero- 
geneity. However, in principle, the inverse is not nec- 
essarily true. Thus, the time for the onset of Fickian 
diffusion is probably only a lower bound for the lifetime 
of dynamic heterogeneities. 

To investigate the onset time we use the trajectories 
generated by an extensive Brownian Dynamics simula- 
tion study of a 80:20 Lennard-Jones binary mixture intro- 
duced by Kob and Andersen ^21]. Briefiy, the potential 
is given by Vap = 46^/3 [(o-a/s/'^)^^ - (cTafs/r)'^], where 
a,(3 e {A,B}, and saa = 1-0, cab = 1-5, esB = 0.5, 
aAA = 1-0, (JAB — 0.8, and asB ~ 0.88 (all the results 
are presented in reduced units where aAA and eAA are 
the units of length and energy, respectively). A total of 
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FIG. 1: The probability of the logarithm of single particle 
displacements P(logj^Q((5r); at T = 0.45 for the A particles. 
The time t is equal to, from left to right, r„g = 160, Tq — 
725, Tnng = 1405, TQ.i = 2990, and tf = 22391 (see text for 
definition of these times). For a comparison we also show, 
as a dotted line, P{logiQ{5r);t) resulting from a Gaussian 
distribution of displacements. 



FIG. 2: Left panel: The probability of the logarithm of single 
particle displacements P(log]^Q((5r); tf) for the A particles for 
T = 1.0, 0.8, 0.6, 0.55, 0.50, 0.47, 0.45, listed from left to 
right. Right panel: the probability distributions from the left 
panel shifted in such a way that {Sr^'^ — 1; for a comparison 
we also show, as a dotted line, P{logj^g{5r);t) resulting from 
a Gaussian distribution of displacements with {Sr^') = 1. 



N = 1000 particles were simulated with a fixed cubic 
box length of 9.4. The details of this study have been 
presented elsewhere . In the present investigation 

we use only some of the temperatures simulated before: 
T = 1.0, 0.9, 0.8, 0.6, 0.55, 0.5, 0.47, and 0.45. The previ- 
ous runs at the temperature T = 0.45 have been extended 
by 60%; the A particles' mean squared displacement at 
the end of the extended runs is about 16. We present the 
results for the A particles only. The results for the B par- 
ticles are qualitatively the same, although the statistics 
is worse due to the smaller number of B particles. The 
temperature dependence is presented by plotting various 
quantities vs. T — Tc where Tc — 0.435 is the crossover 
temperature jlil . l2ll |. This is a convenient way to ex- 
pand the temperature scale and it should not imply an 
endorsement of any particular theoretical approach. 

We start by showing in Fig. the probability distri- 
butions P{\ogiQ{6r)]t) a,tT — 0.45 for the A particles at 
several times characteristic for the relaxation of the sys- 
tem. The first one is the time at which the non-Gaussian 
parameter a2{t) — | (Jr*) / (^Sr^'j — 1 reaches the maxi- 
mum value, r„g. The second one is the a relaxation time, 
Tq,, which is defined in the usual way: is the time at 
which the incoherent intermediate scattering function for 
a wave vector near the peak of the static structure factor 
is equal to 1/e of its initial time value, Fs{k;Ta) = 1/e. 
The third time is the time at which a new non-Gaussian 
parameter 7(i) = ^ (l/5r^) — 1 reaches the 
maximum value, r„„g. We argued in Ref. [19t | that devi- 
ations of P{\ogiQ{Sr);t) from its Fickian shape are most 
evident for times comparable to T„„g. The fourth time 
is the time at which the incoherent intermediate scat- 
tering function for a wave vector near the peak of the 
static structure factor is equal to 10% of its initial time 
value, Fs{k; tq.i) = 0.1. The final time is the onset time. 



Tp, i.e. the time at which the peak of P{\ogiQ{5r);t) 
is equal to the 90% of its value for a Gaussian distribu- 
tion of displacements. For a comparison we also show 
a P{logiQ{Sr);t) resulting from a Gaussian distribution 
of displacements. It is clear from Fig. ^that at shorter 
times, i.e. at r„g, Tq, r„„g, and tq.i the probability dis- 
tributions P{\ogiQ{6r);t) deviate strongly from the shape 
resulting from a Gaussian distribution of displacements. 
While there are still noticeable differences even at Tp , we 
believe that these are small enough to consider Tp the 
onset time for Fickian diffusion. 

In Fig. 121 we show P{\ogiQ{6r);Tp) for the A parti- 
cles for T = 1.0, 0.8, 0.6, 0.55, 0.50, 0.47, and 0.45. It 
should be noted that with decreasing temperature the 
probability distributions at Tp shift toward larger dis- 
placements. In other words, mean squared displacement 
at the onset of Fickian diffusion increases with decreasing 
temperature. The right panel indicates that the shape of 
P{logiQ{Sr);Tp) is temperature-independent and, there- 
fore, the late-time liquids' dynamics are, up to rescaling 
of the time and distance scales, similar. 

Fig. 131 presents our main result: comparison of the 
temperature dependence of of the onset time for Fick- 
ian diffusion, Tp, with that of other characteristic times. 
Wc find that, in the temperature range considered in this 
Letter, the onset time is between 6.5 and 31 r^. More 
importantly, the ratio of the onset time and the a relax- 
ation time grows with decreasing temperature. Interest- 
ingly, the temperature dependence of this ratio becomes 
somewhat weaker with decreasing temperature: it ap- 
pears stronger in the range 0.1 < T — Tc < 1 than in the 
lower temperature range 0.01 < T — Tc < 0.1. 

In Fig. 2|we place the results shown in Figs. [21and|21in 
the context of the time dependence of the mean squared 
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FIG. 3: Temperature dependence of characteristic times 
for the A particles. Triangles: the time at which the non- 
Gaussian parameter 02 (i) reaches the maximum value, r„g. 
Diamonds: the a relaxation time, Tq. Circles: the onset time 
for Fickian diffusion, tf- Inset: temperature dependence of 
the ratio tf/tc. 
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FIG. 4: The time dependence of the mean square displace- 
ment for the A particles for T = 1.0, 0.8, 0.6, 0.55, 0.50, 0.47, 
and 0.45 listed from left to right. The symbols are placed 
at different characteristic times. Squares: the time at which 
the standard non-Gaussian parameter reaches the maximum 
value. Trig. Triangles: the a relaxation time, Tq. Circles: the 
onset time for Fickian diffusion, tf- 



displacement. On the time scale of T„g the mean squared 
displacement has not yet reached the linear dependence 
on time and thus the diffusion is obviously non-Fickian. 
Moreover, on the time scale of Tq, the mean squared dis- 
placement is, at most, at the borderline of the linear time 
dependence. On the other hand, the onset time, Tp, oc- 
curs well within the regime of apparent linear time de- 
pendence of the mean squared displacement. Note that 
there is an important practical message from Fig. ^ if 
one monitors only the time-dependent mean square dis- 
placement, one can significantly underestimate the length 
of the run necessary to achieve Fickian diffusion. 

Having identified the Fickian crossover time, Tp, we 
can define a characteristic length scale, the root mean 
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FIG. 5: Circles (left vertical axis): the root mean squared 
displacement at the onset time, {(^Sr^{TF)'))^^'^ ■ Triangles 
(left vertical axis) : the square root of the product of the self- 
diffusion coefficient and the a relaxation time (multiplied by 
20 for convenience), 20 (Dra)^^^. Squares (right vertical axis): 
the ratio of the onset time and the a relaxation time, tf/tc- 
All quantities are plotted vs. the a relaxation time; all the 
data pertain to the A particles. Dashed lines indicate scaling 

relationships ((^^^(t'f))) 



' ^ oc Ta^^ and tf/tq (x t^'^^. 
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squared displacement at the onset time, ((i5r^(rF))) 
It follows from Figs. [21 and 0| that this length increases 
with decreasing temperature. Our characteristic length 
should be related to the so-called Fickian crossover length 
I* introduced by Berthier et al. The latter length was 
defined through the wavevector dependence of the relax- 
ation time of a supercooled liquid. Roughly speaking, I* 
is the length scale on which diffusion is Fickian on all 
time scales. The prediction of Ref. 23] was that this 
length scale changes with temperature as the square root 
of the product of the self-diffusion coefficient and the a 

relaxation time, I* oc (DTa)^^"^ In Fig- El we com- 

1/2 

pare temperature dependence of (((5r^(rF))) to that 
of of I* (note that we plot these length scales vs. the 
a relaxation time; this is in the spirit of Refs. 0[ 
in which it is argued that the glass transition is a mani- 
festation of a zero-temperature critical point). The root 
mean squared displacement at the onset time grows with 

increasing and at the lowest temperatures there is an 

1/2 

apparent scaling relationship, {{Sr'^ (rp)')) oct^'^^. In 
1/2 

contrast, {Dto) is initially temperature- independent. 

This is due to the fact that the Stokes-Einstein relation is 

violated only for > 1 {i.e. for T < 0.8) [2^. However, 

at longer a relaxation times {i.e. at lower temperatures) 
1/2 

{Dva) has a temperature dependence similar to that 
1/2 

of {{SP{tp))) ' . Finally, we show in Fig. jsl that at 
longer a relaxation times {i.e. at lower temperatures) 
the ratio Tp / Ta appears to grow with increasing a relax- 
ation time as tp/to oc t^-^^. It is not clear whether the 
scaling relations indicated in Fig. [3 have any deeper sig- 
nificance. It could be even argued that if they continue 



4 



for another 7 or 8 orders of magnitude of Tq {i.e. up to 
Ta comparable to that at the laboratory glass transition 
temperature) , the resulting Tp would be greater than the 
longest experimentally observed heterogeneity lifetime. 

To summarize, wc proposed a quantitative definition 
of the onset time for Fickian diffusion and investigated 
its temperature dependence in a Lennard-Jones binary 
mixture. We found that the onset time is considerably 
longer than the a relaxation time and, more importantly, 
it increases faster with decreasing temperature than the 
a relaxation time. Our definition of the onset time re- 
lies upon one particular indicator of Fickian diffusion, 
the probability distribution of the logarithm of single- 
particle displacement, P{logiQ{Sr);t), and upon adopting 
a particular numerical criterion for the onset of Fickian 
diffusion, peak height equal to the 90% of its Fickian 
value. This procedure seems reasonable in that it results 
in non-Fickian motion being present only at tempera- 
tures at and below T « 1.0. This temperature has been 
identified before as so-called onset temperature for slow 
dynamics To test the robustness of our main re- 

sult we also tried using a different indicator of Fickian 
diffusion: the new non-Gaussian parameter that we in- 
troduced recently 0|. In this approach wc defined the 
onset time for Fickian diffusion to be the time at which 
the new non-Gaussian parameter is equal to 1/3. This 
particular numerical value results in non-Fickian motion 
being present only at temperatures at and below T « 0.8 
[isf. The resulting onset times are somewhat shorter 
than the ones presented in this Letter. However, the tem- 
perature dependence of the onset time defined using the 
new non-Gaussian parameter is similar to that of the on- 
set time defined using P(log;^o(^''); More interestingly, 
we found that the shapes of P{logiQ{Sr);t) at the onset 
times defined using the new non-Gaussian parameter are 
very similar and, in particular, the height of the peak 
is approximately temperature independent and equal to 
85% of its value for Fickian diffusion. 

Finally, we would like to point out that the results 
presented here violate the time-temperature superposi- 
tion principle: in order to superimpose the probability 
distributions P{\ogif^{6r);t) shown in the left panel of 

1/2 

Fig. 121 we have to shift logjQ((5r) by log^Q (((5r^(Ti?))) 
The more usual shift procedure, agreeing with the time- 
temperature superposition, would involve the a relax- 
ation time rather than the onset time Tp which has tem- 
perature dependence different from r^. 
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